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Abstract — Heptathela spiders of Okinawajima Is. were taxonomically revised using DNA sequence data. 
An analysis of mt-COI gene sequences revealed two clades (HY and HH), one of which (HY) was identified 
as Heptathela yanbaruensis because HY included specimens from Yona, the type locality of this species. To 
evaluate gene flow between HY and HH, the sequences of the nuclear 28S-rRNA and histone H3 genes were 
analyzed for specimens collected from two sympatric localities. In the two nuclear genes of HY and HH, all 
specimens had 5 homozygous positions with bases that differed between HY and HH. This indicates that 
there has been no gene flow between HY and HH, which is solid evidence for their reproductive isolation. 
Thus, it was concluded that HY and HH should be recognized as independent species, and HH is described 
as a new species under the name Heptathela hellos n. sp. There were no good morphological diagnostic fea¬ 
tures other than the body size of the adult males of the two species {yanbaruensis < hellos)', thus, females 
of these two species can be distinguished only by DNA sequence data. 

Key words — Heptathela hellos, new species, reproductive isolation, species delimitation, species boundary, 
DNA, nuclear 


Introduction 

A total of nine species have been described from Japan in 
the genus Heptathela (Ono 2009) and one of these species, 
Heptathela yanbaruensis Haupt 1983, is found on Okinawa¬ 
jima Is. However, previous studies may have underesti¬ 
mated the true diversity of this genus because spiders with 
limited dispersal capabilities may contain cryptic species 
without distinctive morphological characters. In fact, DNA 
sequence data showed that Aptostichus simus, a trapdoor 
spider with limited dispersal capability, is geographically 
subdivided and divergent in the absence of morphological 
differentiation (Bond et al. 2001). A similar geographical 
divergence was detected in Ryuthela (Tanikawa 2013a, b), 
a sister group of Heptathela. Therefore, we collected many 
specimens of Heptathela yanbaruensis from its entire distri¬ 
bution range and looked for the presence of cryptic species 
using DNA sequence data. 

This study was performed in the following three steps. 
First, we examined partial sequences of the mitochondrial 
cytochrome oxidase subunit I (mt-COI, DNA barcoding 
gene) from all the specimens, with the goal of identifying 
phylogenetic groups and their distribution ranges. Second, 
collected many specimens from localities in which spiders 
of different phylogenetic groups are sympatric and analyzed 


partial sequences of the nuclear 28S-rRNA (28S) and nu¬ 
clear histone H3 (H3) genes to examine potential gene flow. 
Finally, we investigated possible morphological differences 
between the groups. 

Materials and methods 

Sampling. Heptathela specimens were collected from 
the northern part of Okinawajima Is. (Fig. 1), which is the 
distribution range of H. yanbaruensis (Haupt 2003). To 
cover the whole distribution range, collection was con¬ 
ducted in a wider area than the previously defined range. 
Three or more adult spiders per locality were generally col¬ 
lected, but juvenile spiders were also included if adults were 
not available in sufficient numbers. Female and juvenile 
specimens were preserved in 99.5% ethanol at 4°C. Male 
adult spiders were preserved in 75% ethanol at room tem¬ 
perature, but the right fourth leg was dissected and pre¬ 
served in 99.5% ethanol at 4°C for molecular analysis. 
Young males were brought back alive to the laboratory and 
reared until becoming adults. 

DNA extraction, polymerase chain reaction and se¬ 
quencing. Genomic DNA was extracted from muscle of 
legs of large individuals or the whole cephalothorax of 
small individuals using DNeasy Blood & Tissue kit 
(Qiagen, Inc., Germantown, MD) or FavorPrep Tissue 
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Fig. 1. Collection localities of Heptathela yanbaruensis and H. sp. Specimens were identified by mt-COI sequence data. Open circle: H. 
yanbaruensis; filled circle: H. sp.; arrows A and B: sympatric localities. H. sp. is described as H. helios n. sp. in this study. 


Genomic DNA Extraction Mini Kit (Favorgen Biotech 
Corp., Ping-Tung). Mt-COI was amplified using the primer 
combination CBl: 5-TAT GTA CTA CCA TGA GGA 
CAA ATA-3' (Jermiin & Crozier 1994) with HCOI-2198: 
5-TAA ACT TCA GGG TGA CCA AAA AAT CA-3' 
(Folmer et al. 1994). The reactants were initially denatured 
for 2 min at 90°C, proceeded with 40 cycles of 15 sec at 
90°C, 20 sec at 50°C, 4 min at 72°C. 28S was amplified 
using the primer combination ZXl: 5'-ACC CGC TGA 
ATT TAA GCA TAT-3' with AS80P1: 5-AGA GCC 
AAT CCT TGT CCC GA-3' (Bond & Hedin 2006). The 
reactants were initially denatured for 2 min at 94°C, pro¬ 
ceeded with 40 cycles of 30 sec at 94°C, 20 sec at 50°C, 2 
min at 72°C. H3 was amplified using the primer combina¬ 
tion H3aF: 5-ATG GCT CGT ACC AAG CAG ACV GC- 
3' with H3aR: 5-ATA TCC TTR GGC ATR ATR GTG 
AC-3^ (Colgan et al. 1998). The reactants were initially de¬ 
natured for 2 min at 94°C, proceeded with 30 cycles of 30 
sec at 94°C, 20 sec at 50°C, 15 sec at 72°C. PCR products 
were purified using the ExoSAP-IT (GE Healthcare Bio- 
Sciences, Co. Ltd., Buckinghamshire, England). The puri¬ 
fied PCR products were sequenced using the BigDye 
terminator cycle sequencing kit (ver.3.1) using the primer 
HCOI-2198 (mt-COI), ZXl (28S), or H3aF (H3) and ana¬ 
lyzed on ABI 3100 or ABI 3130x1 automated DNA 
sequencer (Applied Biosystems, Foster City, CA). 
Chromatograms were checked by eye using MEGA version 
5.05 (Tamura et al. 2011). Sequence alignments were done 
by MUSCLE (Edgar 2004) in MEGA. Overall mean p- 
distance was calculated using MEGA by averaging the num¬ 
ber of base differences per site over all sequence pairs. The 
phylogenetic tree was constructed by Neighbor Joining 


method using MEGA. 

Morphology. The morphological characters were exam¬ 
ined under stereoscopic microscope, M3Z (Wild Heerbrugg 
AG, Heerbrugg, Switzerland), photographs were taken by 
EOS D60 or EOS Kiss X7 with EFIOO mm macro lens and 
MT-24EX macro twin flash (Canon Inc., Tokyo, Japan) or 
attached to stereoscopic microscope. The relationship be¬ 
tween width of carapace and distance from sympatric locali¬ 
ties were compared for checking the occurrence of 
“character displacemenf’. All measurements were given in 
mm. 

Deposit. Nucleotide sequence data are available in the 
DDBJ/EMBL/GenBank databases. The type specimens 
designated in this study are deposited in the collection of the 
Department of Zoology, National Museum of Nature and 
Science, Tokyo (NSMT). 

Results 

In total, 210 Heptathela specimens from 47 localities 
(Fig. 1, Appendix) were sampled and mt-COI partial se¬ 
quences of 578 bp were obtained. Haplotype names and ac¬ 
cession numbers are shown in the Appendix. In the 
unrooted NJ tree, two clades were detected (Fig. 2A, B). 
Specimens in clade A were identified as Heptathela 
yanbaruensis (HY) because this clade included specimens 
from Yona, the type locality of the species. The remaining 
specimens (those of clade B) were designated as the spiders 
in question (HH). Of the 210 collected specimens, 175 were 
identified as HY and 35 as HH. The uncorrected p- 
distances (number of base differences per site) were 0.002- 
0.130 among HY, 0.002-0.106 among HH, and 0.135-0.196 
between HY and HH. The distribution range of HH was 
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0.02 

Fig. 2. Unrooted bootstrap consensus tree obtained by the Neighbor Joining method for the specimens examined. Scores at two nodes are 
bootstrap values (1000 replicates). Specimens of clade A are identified as Heptathela yanbaruensis and those of clade B as H. helios n. sp. 
in this study. The scale bar shows the substitution per site. 


limited to the northernmost part of Okinawa)ima Is. (Fig. 1). 

A total of 63 specimens were obtained from two 
sympatric localities (Fig. lA, B), of which 58 individuals 
were determined to be HY and 5 to be HH based on the mt- 
COI sequence. In the nuclear 28S partial sequence (627 bp), 
all the 58 HY individuals were homozygous for T at posi¬ 
tions 422 and 603, while all 5 HH individuals were homozy¬ 
gous for C. In the nuclear H3 sequence (303 bp), all 58 HY 
individuals were homozygous for C at positions 103 and 
196, and G at position 235; while all 5 HH individuals were 
homozygous for T at positions 103 and 196, and A at posi¬ 
tion 235. 

No morphological gaps were found in the male palpal 
structure of HY (Fig. 4H-J) and HH (Fig. 4E-G), which is 
the only useful morphological feature for identification of 
Heptathela spiders. The only clear morphological differ¬ 
ence between HY and HH was the body size of adult males: 
the range of the carapace width of adult male HY specimens 
(n= 15) was 3.48-4.25, while that of HH (n=5) was 4.63- 


5.50 (Fig. 3B), though the male HH specimens used in this 
study were only 5. The difference between the two was sig¬ 
nificant (P < 0.005, Mann-Whitney U-test). In female speci¬ 
mens, the average carapace width of HY (3.80, n=61) was 
smaller than that of HH (4.86, n=8), but the ranges over¬ 
lapped to some extent (HY, n = 61: 3.08-4.75, HH, n = 8: 
4.38-5.89) (Fig. 3A). The difference between HY and HH 
were significant again (P< 0.001, Mann-Whitney U-test). 

Discussion 

Two nuclear genes of HY and HH had 5 positions where 
all the specimens exhibited homozygosity, but the bases of 
HY and HH were different. If intercross between HY and 
HH had occurred, these positions should have been het¬ 
erozygous. Therefore, there has been no gene flow between 
HY and HH, which is strong evidence for their reproductive 
isolation. Consequently, it is concluded that HY and HH 
are independent species and HH is described as a new spe¬ 
cies in the following section. 
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Distance from sympatric populations (km) Distance from sympatric populations (km) 

Fig. 3. Relationship between the earapaee width (mm) and the distanee from sympatrie populations (km); + means northeast direetion 
and — means southwest direetion. A: female, B: male. Open eireles: HY {Heptahela yanbaruensis), solid eireles: HH (deseribed as H. 
helios n. sp in this srudy). 


A previous study (Haupt 2003) revealed that females of 
Heptathela or related spiders attraet males with a kind of 
pheromone and that males produee vibrational signals using 
their palps when approaehing female burrows. Therefore, 
reproduetive isolation between HY and HH may be retained 
by the differenee in ehemieal substanees made by females 
and/or the differenee in male courtship signals. The mating 
season of these two species may also differ because mature 
HY males were found in late October to early December, 
while no adult HH males were found in this season. All the 
mature male specimens of HH used in this study were col¬ 
lected as juveniles and reared until becoming adults, there¬ 
fore the maturation period of HH males is still unknown. 

The large p-distance for the mt-COI sequences between 
HY and HH suggests a long history after speciation of these 
two species. The distribution range of HH is currently re¬ 
stricted to the northernmost part of Okinawajima Is., but this 
range may have been wider in the past, and probably north¬ 
ward, when the land area was larger than present. 

Distributional map of HH and HY indicates that the 
distributional zone of overlap between the two is rather nar¬ 
row (Fig. 1). It is likely that selective forces in evolution of 
some morphological characters is different between 
sympatric populations where HH and HY co-occur and 
allopatric populations where only one of the two forms oc¬ 
curs if there is some interspecific competition between HH 
and HY. In this sense, it is interesting that female body size 
of HY seems to be decreased in the sympatric zones, while 
that of HH appears to be increased in the same zone than 
nearby non-sympatric populations of respective species 
(Fig. 3A). Thus, it seems that this situation in the zone of 
overlap indicates occurrence of “character displacement” in 
female body size of the two species, although, unfortu¬ 
nately, we have no solid answer whether this phenomenon 
is observed also in males since no adult males of HH have 
been obtained from the sympatric populations (Fig. 3B). 


Description of New Species 
Heptathela helios n. sp. 

[Japanese name: Kunjan-kimuragumo] 

(Fig. 4A-G) 

Type series. All the type specimens were collected from 
Kunigami-son, Okinawajima Is., Japan, by Akio Tanikawa 
unless noted otherwise. Holotype: Ic/', Oku (N26.84274 
E128.28737), 26-V-2010 (collected as juvenile and became 
adult on 9-IX-2012) (NSMT-Ar 12851). Paratypes: 2c^, 
same data as holotype (became adults on 15-X-2010 and 5- 
XII-2012) (NSMT-Ar 12852-12853); 1^, Oku (N26.83633 
E128.27544), 30-V-2007 (NSMT-Ar 12854); 2c^, Oku 
(N26.83436 E128.27615), 26-V-2010 (collected as juveniles 
and became adults on 20-X-2010 and 13-1-2013) (MSNT- 
Ar 12855-12856); 1^, Oku (N26.83436 E128.30238) 
(NSMT-Ar 12857), 2^, Benoki (N26.79985 E128.26064), 
lO-XII-2012 (NSMT-Ar 12858-12859); 1^, Benoki 

(N26.79797 E128.28802), ll-XII-2012, Booppa Petcharad 
leg. (NSMT-Ar 12860); 1^, Ada (N26.75557 E128.30072), 
13-VI-2006 (NSMT-Ar 12861), 1^, same locality, 11-XII- 
2012 (NSMT-Ar 12862). 

Etymology. Arbitrary naming. Helios is the personifica¬ 
tion of the sun in Greek mythology. “Kunjan” of the 
Japanese name is an alias name of “Kunigami” which refers 
to the northern part of Okinawajima Island. 

Diagnosis. This species can be distinguished from H. 
yanbaruensis by the following unique combinations of mt- 
DNA COI gene nucleotide site substitutions; C (21), C 
(262), C (266), A (384), C (558). The body size of adult 
males of H. helios is larger than that of H. yanbaruensis, but 
the male palps of these two species are quite similar (com¬ 
pare Fig. 4E-G and 4H-J), which cannot be used for identi¬ 
fication. Females of H. helios also cannot be distinguished 
from H. yanbaruensis by morphological features. 

Description. Based on the holotype ^ (NSMT-Ar 
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Fig. 4. Heptathela helios n. sp. (A-G) and H. yanbaruensis (H-J). A, habitus of female (paratype from Benoki, NSMT-Ar 12858); B, 
habitus of male (holotype from Oku, NSMT-Ar 12851); C, female genital area, ventral view (paratype from Benoki, NSMT-Ar 12858); D, 
spermatheeae, dorsal view (paratype from Ada, NSMT-Ar 12861); E-G, male palp, distal view (E: holotype from Oku, NSMT-Ar 12851; 
F-G: pratypes from Oku, NSMT-Ar 12852, 12855; H-J: from Nago). Seales = 5 mm (A-B), 1 mm (C-F). 
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12851, Fig. 3B) and one of female paratypes (NSMT-Ar 
12858, Fig. 3A). Coloration and markings in alcohol. 
Female and male. Carapace pale brown, mottled with dark 
brown, head region anteriorly darker. Abdomen pale brown, 
tergites dark brown. 

Measurements. J'/^, measurements in parentheses indi¬ 
cate the range among type series. Body 12.77 (10.15- 
13.38)/14.62 (10.31-14.62) long. Carapace 5.91 (5.13- 
6.13)/6.50 (4.69-6.50) long; 5.13 (4.63-5.50)75.80 (4.38- 
5.80) wide. Length of legs [tarsus + metatarsus + tibia + 
patella + femur = total]: I, 2.63 + 5.00 + 3.84 + 2.38 

+ 5.38 = 19.20 / 1.69 + 3.00 + 2.81 + 2.63 + 4.91 = 

15.04; II, 3.00 + 5.31 + 3.69 + 2.38 + 5.19 = 19.57/ 

1.88 + 3.13 + 2.56 + 2.63 + 4.69 + 14.89; III, 3.50 + 

6.25 + 3.50 +2.44 + 5.00 = 20.69/2.13 + 3.63 + 2.50 
+ 2.63 + 4.44 = 15.33; IV, 4.19 + 8.19 + 5.00 + 2.66 

+ 6.69 = 26.73 / 3.06 + 5.88 + 3.88 + 3.00 + 6.31 = 

22.13. Abdomen 5.13 (4.25-6.06)78.25 (4.81-8.25) long, 

3.88 (3.13-4.69)76.10 (3.25-6.30) wide. 

Carapace longer than wide [length 7 width 1.15 (1.08- 
1.15)71.12 (1.07-1.22)]. Eye area wider than long [length/- 
width 0.85 (0.83-0.88)70.79 (0.79-0.90)]. Fang furrow of 
chelicera with 12 (7-12)713 (11-13) anterior marginal teeth 
and no retromarginal teeth. Labium wider than long 
[length/width 0.48 (0.38-0.50)70.53 (0.48-0.53)]. Sternum 
longer than wide [length/width 2.15 (1.78-2.18)71.92 (1.65- 
2.00)]. Length of leg I/length of carapace 3.25 (3.16-3.39)/ 
2.31 (2.07-2.34). Male palp as in Fig. 4E-G. Abdomen 
longer than wide [length/width 1.32 (1.29-1.36)71.35 (1.27- 
1.48)], female genital area and spermatheca as in Fig. 4C- 
D. 

Distribution. Northernmost part of Okinawajima Is., 
Nansei Isis., Japan. 
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Appendix. Localities, haplotype or genotype names, and DDBJ/EMBL/GenBank accession numbers of sequence data analyzed in the study. 


Heptathela sp. is described as Heptathela helios n. sp. 

in this paper. A and B indicate the sympatric localities corresponding to Fig. 1. 

Locality 

N 

COI 

Haplotype 

Accessin No. 

28S Genotype 

Accession no. 

H3 Genotype 

Accession no. 

Heptathela sp. 

N26.85868 E128.26139 

1 

OKHH05 

AB918042 

— 

— 

— 

— 

N26.85868 E128.26139 

1 

OKHH09 

AB918045 

— 

— 

— 

— 

N26.84274 E128.28737 

1 

OKHH13 

AB918049 

— 

— 

— 

— 

N26.84120 E128.28388 

1 

OKHH26 

AB918062 

— 

— 

— 

— 

N26.83839 E128.25859 

1 

OKHH04 

AB918041 

— 

— 

— 

— 

N26.83633 E128.27544 

1 

OKHHIO 

AB918046 

— 

— 

— 

— 

N26.83532 E128.28658 

2 

OKHH06 

AB918043 

— 

— 

— 

— 

N26.83532 E128.28658 

1 

OKHH07 

AB918044 

— 

— 

— 

— 

N26.83436 E128.30238 

2 

OKHH02 

AB918039 

— 

— 

— 

— 

N26.83436 E128.27615 

1 

OKHH20 

AB918056 

— 

— 

— 

— 

N26.81580 E128.30493 

5 

OKHH15 

AB918051 

— 

— 

— 

— 

N26.81580 E128.30493 

1 

OKHH23 

AB918059 

— 

— 

— 

— 

N26.81286 E128.31701 

1 

OKHH03 

AB918040 

— 

— 

— 

— 

N26.81286 E128.31701 

1 

OKHHll 

AB918047 

— 

— 

— 

— 

N26.81254 E128.29852 

1 

OKHH17 

AB918053 

— 

— 

— 

— 

N26.80313 E128.29236 

1 

OKHHOl 

AB918038 

— 

— 

— 

— 

N26.79985 E128.26064 

3 

OKHH14 

AB918050 

— 

— 

— 

— 

N26.79985 E128.26064 

1 

OKHH18 

AB918054 

— 

— 

— 

— 

N26.79985 E128.26064 

2 

OKHH19 

AB918055 

— 

— 

— 

— 

N26.79797 E128.28802 

2 

OKHH16 

AB918052 

— 

— 

— 

— 

A: N26.77713 E128.26423 

1 

OKHH25 

AB918061 

OKHH2801 

AB918102 

HH02 

AB918106 

B: N26.75557 E128.30072 

1 

OKHH12 

AB918048 

OKHH2801 

AB918102 

HHOl 

AB918105 

B: N26.75557 E128.30072 

1 

OKHH21 

AB918057 

OKHH2801 

AB918102 

HH04 

AB918108 

B: N26.75557 E128.30072 

1 

OKHH22 

AB918058 

OKHH2801 

AB918102 

HH03 

AB918107 

B: N26.75557 E128.30072 

1 

OKHH24 

AB918060 

OKHH2801 

AB918102 

HH04 

AB918108 

Heptathela yanbaruensis 

N26.83633 E128.27544 

1 

OKHY203 

AB778252 

— 

— 

— 

— 

N26.81074 E128.23613 

4 

OKHY114 

AB918075 

— 

— 

— 

— 

N26.79795 E128.25816 

9 

OKHY114 

AB918075 

— 

— 

— 

— 

N26.79303 E128.31055 

4 

OKHY104 

AB918066 

— 

— 

— 

— 

N26.78692 E128.22362 

2 

OKHY202 

AB918097 

— 

— 

— 

— 

A: N26.77713 E128.26423 

1 

OKHY203 

AB778252 

OKHY2809 

AB918103 

HY02 

AB918110 

A: N26.77713 E128.26423 

1 

OKHY203 

AB778252 

OKHY2801 

AB778049 

HY06 

AB918114 

A: N26.77713 E128.26423 

1 

OKHY203 

AB778252 

OKHY2801 

AB778049 

HY07 

AB918115 

A: N26.77713 E128.26423 

1 

OKHY203 

AB778252 
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